The Inflaton Field as Self-Interacting Dark Matter in the Braneworld Scenario 
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A unified model is developed within the context of the braneworld paradigm, where a single scalar 
field can act as both the inflaton field in the very early universe and also as strong, self-interacting 
dark matter in the post-inflationary universe. Reheating proceeds due to the overproduction and 
subsequent evaporation of primordial black holes. Observational constraints, most notably from 
gravitational waves, are satisfied if the probability of PBH formation is sufficiently high. 
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The recent advances in astronomical observations are 
casting new light on old problems in cosmology. Two 
of the most fundamental questions of interest today are: 
(i) the origin of the scalar inflaton field responsible for 
the accelerated, inflationary expansion of the very early 
universe and: (ii) the nature of dark matter in the uni- 
verse. Both questions have important implications for 
our understanding of large-scale structure formation . 

New ideas about inflationary cosmology are now 
emerging with the development of the braneworld sce- 
nario, where our observable universe is viewed as a do- 
main wall embedded within a higher-dimensional space 
|^|,D . A striking feature of this scenario is the presence of 
a quadratic density term in the Friedmann equation [Q. 
Under quite general conditions, this term allows steep 
scalar field potentials to support an inflationary epoch 
that would otherwise be impossible in standard cosmol- 

ogy||@-@. 

An attractive feature of steep inflationary models is 
that the universe can be naturally reheated by the pro- 
cess of gravitational particle production, where particles 
are produced quantum mechanically at the end of infla- 
tion due to the time-variation of the gravitational field 
Jl0| . ff the inflaton is stable and is able to survive through 
to the present epoch, it may represent a possible candi- 
date for the cosmological constant or quintessence field 
lyUHI that has been proposed to account for the high 
redshift type fa supernova data flif| , |i"2[ . 

On the other hand, Barshay and Kreyerhoff first pro- 
posed that a metastable inflaton field could be identified 
with the present-day cold dark matter (CDM) |f3|. By 
employing renormalization-group techniques to calculate 
radiative corrections to a scalar field potential involving 
a quartic self-interaction these authors found that 
inflation could proceed in the vicinity of the potential's 
maximum and would end as the field reached a well- 
defined minimum. Fluctuations about this minimum re- 
sult in massive inflaton quanta with an energy density 
that is sufficiently high that they can provide a signifi- 
cant fraction of the CDM at the present epoch. 



In this paper, we consider whether the inflaton field 
can be identified with the cold dark matter in the uni- 
verse fl3|| , within the specific context of the braneworld 
scenario and the so-called 'strong, self-interacting scalar 
field dark matter' (SFDM) hypothesis that has recently 
been developed by two of the authors Jf5|-pf[ (see also 
p2| ) . The key idea of the SFDM scenario is that the dark 
matter responsible for structure formation in the universe 
is a real scalar field, <&, that is minimally coupled to Ein- 
stein gravity and has self-interactions parametrized by a 
potential energy of the form 



V{9) = V [cosh(aV^) - 1] , 



(f) 



where Vq and a are the only two free parameters of the 
model, k = 8irG and we employ natural units such that 
Ti = c = I. The effective mass of the scalar field is given 
by m% — kqVqo 2 . A minimal coupling to gravity avoids 
the strong restrictions imposed by the equivalence prin- 
ciple on scales of the order of the solar system. 

The advantage of the SFDM model is that it is insen- 
sitive to initial conditions and the field behaves as CDM 
once it begins to oscillate around the minimum of its po- 
tential. In this case, it can be shown Jl^,|f^] that the 
SFDM model is able to reproduce all the successes of 
the standard ACDM model above galactic scales. Fur- 
thermore, it predicts a sharp cut-off in the mass power 
spectrum due to its quadratic nature, thus explaining the 
observed dearth of dwarf galaxies, in contrast with the 
excess predicted by high resolution N-body simulations 
with standard CDM Jf^,|2^]. The strong self-interaction 
of the scalar field results in the formation of solitonic ob- 
jects called 'oscillatons', which have a mass of the order 
of a galaxy but do not exhibit the cusp density profiles 
characteristic of standard CDM [pjj|~pf|. The best-fit 
model to the cosmological data can be deduced from the 
current densities of dark matter and radiation in the uni- 
verse and from the cut-off in the mass power spectrum 
that constrains the number of dwarf galaxies in clusters. 
The favoured values for the two free parameters of the 
potential (|l|) are found to be |Tq] : 



f 



a ~ 20.28, 
Vb ~ (3 x IQ- 27 M 4 



(2) 
(3) 



where M 4 = G x l 2 w 10 5 gis the four-dimensional 
Planck mass. This implies that the effective mass of the 
scalar field should be ~ 9.1 x 10~ 52 M 4 = 1.1 x 10~ 23 
eV. 

An important feature of the potential (|l]) is that it 
is renormalizable and exactly quantizable, although it 
is presently unknown whether it originates from a fun- 
damental quantum field theory |p4j . Furthermore, the 
scattering cross section by mass of the scalar particles, 
<72^2 /to*) can be constrained from numerical simulations 
of self-interacting dark matter that avoid high-density 
dark matter halos p5| . This effectively constrains the 
renormalization scale, A$, of the potential to be of the 
order of the Planck mass, A* ~ 1.93M 4 = 2.15 x 10 19 
GeV HU . Such a value is indicative of a possible funda- 
mental origin for the scalar field, which in turn suggests 
that the strongly, self-interacting scalar field dark mat- 
ter may also have been present during the inflationary 
epoch p8| . However, we do not perform quantum field 
theoretic or semi-classical calculations in this paper. 

In view of the high energy scales associated with the 
very early universe, it is natural to assume that the 
scalar field was initially displaced from its global min- 
imum, n a 2 & 2 ^> 1. In this limit, the potential is well 
approximated by an exponential function, but from Eq. 
(H), the self-coupling is too large to support inflation- 
ary expansion in a standard cosmological setting. On 
the other hand, such a potential can drive inflation suc- 
cessfully within the braneworld scenario |7|||. This fol- 
lows because the Friedmann equation is modified due to 
the motion of our observable universe (the domain wall) 
through the higher-dimensional 'bulk' spacetime. In par- 
ticular, in the second Randall-Sundrum scenario || , the 
Friedmann equation is given by 



(4) 



when appropriate conditions are satisfied, where H = 
a/a is the Hubble parameter, a is the scale factor of the 
universe, p is the energy-density of the inflaton field (as- 
sumed to be confined to the brane) and Af, is the tension 
of the brane. Conventional Einstein gravity is recovered 
in four dimensions when the energy density is signifi- 
cantly lower than the brane tension, p <C A;,. However, at 
high energies, the quadratic correction implies that the 
expansion rate of the brane is enhanced relative to what 
it would be in a universe governed by Einstein gravity || . 
Thus, the friction acting on the scalar field is increased 
and inflation driven by a potential of the form ([!]) is then 
possible at sufficiently high energies even though a 2 1. 

Braneworld inflation driven by such a potential has 
been studied in Refs. |7|||. Recalling the main results, 
the COBE normalization of the cosmic microwave back- 
ground (CMB) power spectrum J2(| relates the value of 



the brane tension to the scalar field self-coupling such 
that A^ 4 a 3 / 2 w 10 15 GeV. For the favoured value of the 
latter, as implied by Eq. (||), we deduce that 

A b ~ (6 x 10~ 7 M 4 ) 4 = 2.88 x 10 51 GeV 4 . (5) 

For these given values of {a, A;,}, the magnitude of 
the potential energy at the end of inflation is V en d — 
(3.2 x 10~ 6 M 4 ) 4 = 2.33 x 10 54 GeV 4 and this implies 
that ^cnd ~ 2M 4 , thus justifying the exponential ap- 
proximation to the potential (|lj) during the inflationary 
era. In comparison with the canonical V = A<I> poten- 
tial, where COBE normalization implies the dimension- 
less parameter A as 10~ 13 should be very small, our model 
requires the small dimensionful number for Vo given in 
Eq. (||), although this constraint does not arise from the 
CMB. Furthermore, the magnitude of the potential N e- 
foldings before the end of inflation is Vn as (N + l)V en< i, 
implying that the initial value of the scalar field is closer 
to $end than in the quartic model. 

Given the COBE normalization [^6j , the spectral index 
of the scalar fluctuation spectrum is determined to be Q] 



n 



1 



N + l 



0.94, 



(6) 



where N as 70 is the number of e-foldings that elapse 
between the epoch that a given, observable mode crosses 
the Hubble radius during inflation and the end of the 
inflationary epoch. Remarkably, the tilt of the scalar 
perturbation spectrum in this scenario is uniquely deter- 
mined by the number of e-foldings and is independent 
of the parameters in the potential (Q). A spectrum with 
a tilt of this magnitude away from scale-invariance is 
presently favoured by analyses of the CMB power spec- 
trum p7[ . Furthermore, the amplitude of the primordial 
gravitational wave spectrum, At, relative to that of the 
density perturbations, As, can be estimated as H 



A 2 

r = 47T— %■ 
A S 



0.4 



(7) 



implying after COBE normalization that A T as 1.7 x 
10~ 10 . This ratio is also independent of the model's pa- 
rameters and is within the projected sensitivity of the 
Planck satellite. It provides a potentially powerful test 
of the model. 

Inflation ends when the quadratic corrections to the 
Friedmann equation (0) become negligible. Due to the 
steep nature of its potential, the inflaton then behaves 
as a massless field, where its energy density redshifts as 
p$> oc a~ 6 . This is important, because after the ten- 
sor modes have re-entered the Hubble radius, the evo- 
lution of the spectral gravity wave energy density, p g , is 
sensitive to the effective equation of state in the post- 
inflationary universe^] |^,|8[ . It is enhanced (reduced) on 



"The equation of state is assumed to be of the form p — up, 



2 



shorter scales if u > 1/3 (u < 1/3). In general, the 
bound on the gravitational waves imposed by success- 
ful nucleosynthesis, p g < 0.2p ra d, must not be violated 
and this results in an upper limit on the duration of the 
kinetic energy dominated phase H. 

Sahni, Sami and Souradeep |J assume that the evolu- 
tion of the short wavelenth gravitational waves is similar 
to that of conventional cosmology. They then conclude 
that the gravitational wave energy density begins to dom- 
inate the scalar field when the universe has expanded by a 
factor of « 10 5 and, consequently, the universe must 
become radiation dominated before the temperature has 
fallen by this factor. Unfortunately, for the model under 
consideration, the thcrmalizcd temperature of the radia- 
tion produced from gravitational particle production at 
the end of inflation is given by T on( j — 2 x 10 9 GeV, 
whereas the temperature at the epoch when this radia- 
tion dominates the scalar field is T cq ~ (1 — 2) GeV [||, 
corresponding to a redshift of 10 9 . 

An alternative mechanism for reheating the universe 
is therefore required that reduces the duration of the ki- 
netic energy dominated phase. One possibility is through 
the overproduction of primordial black holes (PBHs) that 
subsequently decay into relativistic particles via Hawk- 
ing evaporation 29 - 3|] . In the above inflationary model, 
a fraction, /3o, of the energy density of the universe col- 
lapses into PBHs due to the density fluctuations that 
re-enter the Hubble radius immediately after the uni- 
verse has ceased to accelerate. The PBHs form with 
a mass of the order of the horizon mass at this time 
p2[ and this is given by Af p bh ~ MlH~^ A , where the 
Hubble radius at the end of inflation, H~^ d , is estimated 
from the Friedmann equation (Q) under the assumption 
that p on d ~ Vend ~ 2a 2 Af,. Eq. (|J) then implies that 
-Vfpbh ~ 10 9 M4 and the lifetime of PBHs with this mass is 
icvap ~ {M phh /M 4 ) 3 t P w 10" 16 s, where t P is the Planck 
time. This is sufficiently short that constraints on PBH 
evaporations from primordial nucleosynthesis are satis- 
fied Q33). 

Once formed, the PBHs behave as a pressureless fluid 
and their energy density redshifts as p p hh oc a~ 3 . Thus, 
they are able to dominate the (massless) scalar field be- 
fore they evaporate if 



0o 



l-0o 



> 



pbh 



M 4 



10 



-18 



(8) 



and this change in the effective equation of state occurs 
after the universe has expanded by a factor 



0o 



l-0o 



"end 
Q-dom 



(9) 



where to is a constant barotropic index. For a massless scalar 
field, uj = 1. 



where a^om denotes the scale factor at the onset of PBH 
domination. It follows, therefore, that the PBHs domi- 
nate the scalar field before the backreaction of the gravi- 
tational waves becomes significant if the initial mass frac- 
tion satisfies (3o > 10~ 15 . This is consistent with Eq. (j|) 
and the constraint arising from the integrated gravita- 
tional wave energy density is therefore alleviated since 
the PBH equation of state, to — 0, is 'softer' than that of 
radiation. 

A stronger, and more reliable, constraint on the ini- 
tial PBH fraction can be imposed by requiring that the 
PBHs dominate the universe before the scalar field has 
reached the minimum of its potential and begun to os- 
cillate. During the kinetic dominated regime, the scalar 
field varies as 



$ = $cnd 



A M4 in / a 

47T \ a e „d 



(10) 



and from the estimate of ^end given above, ^end ~ 2M 4 , 
the field reaches the minimum of its potential after the 
universe has redshifted by a factor a/a en d ~ 60. Thus, 
from Eq. (^), the PBHs dominate the universe before 
this point is reached if (3o > 5 x 10~ 6 . 

If the PBHs come to dominate sufficiently early, the 
displacement of the scalar field away from its minimum 
is such that the potential is still well approximated by an 
exponential form at this time. This is the case for $ > 
0.01Af4. Moreover, the standard Friedmann equation is 
valid for <f> < I.8M4, where the quadratic corrections in 
Eq. (|^) become negligible. Since its potential is steep, 
the field subsequently tracks the PBH (fluid) component 
in this regime as in the standard cosmology f3~ij] , where 
its potential and kinetic energies scale at the same rate 
as that of the PBH energy density. More specifically, 
ft$ = 3(l+w)/a 2 and $ 2 /U = 2(1 + uj)/{\ - uj) and this 
implies that the variation of the scalar field with respect 
to the scale factor during tracking is given by [B4[ 



_ _ 3(l+a;) . / a 

* = $to f= ln 

jK a \a t o 



(11) 



where a subscript 't0' denotes the onset of the tracking 
behaviour. This is important because the large value of 
the field's self-coupling, a « 20, implies that the uni- 
verse expands by many orders of magnitude before the 
field reaches its global minimum. In particular, if the 
PBHs dominate the cosmic dynamics for the majority of 
their lifetime, the universe can expand by up to a factor 
of 10 12 before the PBHs evaporate. The rolling of the 
scalar field down its potential during this epoch is only 
A$ « —0.84. The subsequent transition to a radiation 
dominated universe has a negligible effect on the track- 
ing behaviour of the field, and for a wide range of initial 
conditions, the field does not reach its minimum until 
well after the primordial nucleosynthesis era has passed. 
In this case, the nucleosynthesis bounds are not violated 
since a > 5 |35j. In principle, therefore, our model does 



3 



not exhibit the problem of overshooting nor undershoot- 
ing for the initial conditions [fjEf . The post-inflationary 
universe after the nucleosynthesis era would then corre- 
spond to the universe considered in Refs. [p^|-^8|, where 
it was shown that the scalar field can subsequently act 
as dark matter in the universe. 

To summarize, we have found that a scalar field with a 
potential of the form ([!]) can drive an epoch of inflation- 
ary expansion in the braneworld scenario and may also 
act as a candidate for the dark matter at the present 
epoch. The reheating of the universe proceeds via PBH 
domination and evaporation. The model is able to ex- 
plain a variety of astrophysical observations over a wide 
range of scales and, indeed, contains only four, free pa- 
rameters: the effective mass of the scalar field, m$, its 
self-coupling, a, the tension of the brane, and the 
initial mass fraction of PBHs, [3q. It is the insensitivity 
of Eqs. (^J) and (0) to the potential's parameters that al- 
lows us to simultaneously satisfy the constraints arising 
from the mass power spectrum without the need for fine 
tuning. Once fixed by observational constraints, the pa- 
rameters of the potential can remain unaltered through- 
out the history of the universe. 

In this paper, we have considered the potential (Q) 
from a phenomenological perspective with the primary 
aim of determining the region of parameter space con- 
sistent with astrophysical observation. The parameters 
cover a wide range of scales and one question that im- 
mediately arises is whether the form of the potential (|l|) 
and, in particular, the best-fit values for its parameters, 
can be more fully understood from a field-theoretic per- 
spective. This is beyond the scope of the present paper. 
Nevertheless, we have developed a model based on a sin- 
gle scalar field that covers the history of the universe from 
the inflationary era through to the present epoch and it 
is therefore to be expected that the parameters should 
cover a similar range of energy scales. The success of 
the model in simultaneously satisfying a number of fun- 
damental observational constraints provides strong mo- 
tivation for considering its field-theoretic origin in more 
detail. 

We have viewed the initial mass fraction of PBHs as a 
free parameter in the above analysis, but its magnitude 
is determined by the density perturbations at the end 
of inflation. A more detailed calculation of the fluctua- 
tion spectra is therefore required if further insight is to 
be gained, but this involves an extension of the slow-roll 
analyses employed thus far Q and is beyond the scope 
of the present paper. Since the post-inflationary uni- 
verse is initially dominated by a scalar field, the question 
of PBH formation in this scenario is closely related to 
the problem of scalar field collapse and this topic has 
attracted considerable attention in recent years. (For a 
recent review, see, e.g., Ref. |36)). Moreover, in esti- 
mating the limits on the probability of PBH formation, 
we have assumed that standard four-dimensional results 
derived within the context of Einstein gravity are valid. 
This is reasonable since the PBHs form once inflation has 



ended and this occurs when the brane corrections to the 
Friedmann equation ^ have become negligible. 

A potential problem with reheating a braneworld infla- 
tionary universe via PBH evaporations is that the decay 
products may radiate primarily off the brane and into the 
higher-dimensional bulk, thereby rendering the brane ef- 
fectively cold and empty. This remains an open question 
in the second Randall-Sundrum scenario, but it has been 
shown in a related model that most of the Hawking radi- 
ation consists of standard model particles that are indeed 
confined to the brane p7| . 

Finally, there is the question of whether the PBHs leave 
behind stable, Planck-sized relics at the endpoint of their 
evaporation |50 38 1. Although such a possibility is now 
considered unlikely, it is worth remarking that when the 
PBHs dominate the universe before they evaporate, their 
relics can not overdose the universe if M p bh < 10 6 g |3(i| l 
and this bound is satisfied for the above scenario. 

In conclusion, we have proposed a simple, unified 
model of the inflaton and dark matter particles, where 
the same scalar field provides the origin for the primor- 
dial spectrum of density perturbations produced quan- 
tum mechanically during inflation and also plays a cen- 
tral role in forming the cosmological structures that we 
observe today. 
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